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Abstract. We combine ASCA and RXTE data of V1062 Tau to confirm the presence of a 62- min X-ray p.dsaUon 
We show that the pulsation is caused largely by the variation of dense part.al absorption, in keeping with cur 
models of accretion onto magnetic white dwarfs. Further parametrisat.on of the spin pulse is, however, hampered 
by ambiguities in the models. 
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1. Introduction 

The intermediate polar (IP) class of cataclysmic variable 
is characterised by prominent X-ray pulsations resulting 
from accretion onto magnetic white dwarfs (see Patterson 
1994 for a review). Several studies have attempted to use 
spectral analysis of X-ray data to deduce the cause of the 
pulsations (e.g. Rosen etal. 1998; Hellier etal. 1991, de 
Martino etal. 2001). In particular, some studies have ar- 
gued that the pulses are caused by accretion regions mov- 
ing over the white-dwarf limb for parts of the white-dwarf 
spin cycle (e.g. Singh & Swank 1993; Allan etal. 1998), 
whereas other studies have suggested that the prime cause 
is absorption occurring when parts of the accretion flow 
cross the line of sight to the accretion regions (e.g. Rosen 
et al. 1998; Hellier et al. 1991). In general, both effects may 
be important. 

In probing further, we have found that spectral anal- 
ysis often gives ambiguous results. Phase-varying absorp- 
tion must be partial, and often of several different densi- 
ties, while a disappearance over the limb must take into ac- 
count the temperature-stratification of the accretion col- 
umn. The lowest, coolest regions disappear first, which 
can mimic the spectral signature of partial absorption. In 
at least one star, V1025 Cen, we found that ASCA data 
alone were not adequate to distinguish between the dif- 
ferent complex models (Hellier et al. 1998). We have thus 
obtained simultaneous ASCA and RXTE data of an IP, 
aiming to combine the spectral resolution of ASCA with 
the harder response and higher throughput of RXTE to 
reduce the ambiguities. As a target we chose V1062 Tau, 
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with the additional aim of confirming its membership of 
the IP class, given that the X-ray spin pulse was previously 
known from only 6 cycles of data (Remillard etal. 1994). 
With both a spin period (62 min) and an orbital period 
(10 h) that are unusually long, V1062 Tau marks one ex- 
treme of the P S p ln -/ > orb plane (e.g. Hellier 1999), though 
it fits the P s p m * O.lPorb rule-of-thumb that many IPs 
obey. 

2. The data 

We observed V1062 Tau between 1998 February 16 06:38 
UT and 1998 February 18 13:25 UT in a simultaneous ob- 
servation using the RXTE and ASCA satellites (see Bradt 
etal. 1993 and Tanaka etal. 1994 respectively). ASCA 
recorded 60 ks of good data during this interval in the 
energv range 0.6— 10 keV. The count rate was 0.3 s in 
each of the 4 instruments (SIS0, SIS 1 , GIS2, GIS3). RXTE 
recorded 51 ks of data with the PGA instrument with all 
5 PCUs working. To maximise the signal-to-noise we used 
only the top xenon layer of PCA data, which after back- 
ground subtraction using PCABACKEST v2.0c gave an av- 
erage of 19 counts s' 1 in the energy range 2-15 keV. 

In Fig. 1 we show the combined 1-15 keV RXTE and 
ASCA lightcurve of V1062 Tau. This results from nor- 
malizing the data from the different instruments to the 
same count rate and then merging the lightcurves. Since 
the instruments have different spectral responses this is a 
dubious procedure, but for period searching the overrid- 
ing consideration is to maximise the data coverage. Fig. 2 
shows the Fourier transform of the lightcurve. There is a 
clear periodicity at 3704 ± 8 s, consistent with the spin 
period detected by Remillard et al. (1994), thus confirming 



Fig. 1. Part of the combined ASCA and RXTE 1-15-keV lightcurve of V’1062 Tau, in 128-s bins 



that V1062 Tau is an IP. There is also power at the first 
harmonic of the spin period (2w). The roughly sinusoidal 
pulse profile is shown in Fig. 3. 

There is no significant power at the 10-hr orbital period 
reported by Remillard etal. (1994). nor any at the beat 
period between the spin and orbital cycles. We can thus 
assume that V1062 Tau is a purely disc-fed accretor (see, 
e.g., Hellier 1991), which is as expected given the long 
orbital period. 


3. Spectral analysis 

For phase-resolved analysis we divided the data into four 
bins, with phases 0.867-0.180 (maximum); 0.180-0.270 
(fall); 0.270-0.695 (minimum); and 0.695-0.867 (rise) VVe 
then added the ASCA SIS0 and SISl spectra, and also 
the GIS2 and GIS3 spectra, so obtaining three spec- 
tra per phase bin (SIS, GIS and PCA) for simultaneous 
analysis in XSPEC. VVe encountered two problems. First, 
the spectrum of the source showed secular changes after 
1998 February 18 02:00 UT. The source first became 30% 
brighter, then, four hours after the brightening, the spec- 
trum became more absorbed, with the hardness ratio f(6- 
10)/ (2—6) keV in the ASCA data] increasing from 0.25 to 



0.50. Thus, for phase- resolved analysis, we excluded the 
data after 02:00 UT. 





Secondly, we found discrepancies in the normalisations 
between the SIS, GIS and PGA spectra, with th5 PCA 
spectra being systematically ~ 20% brighter than the GIS 
spectra, which were ^ 10% brighter than the SIS spectra. 
Some part of this is accounted for by the different sam- 
pling of the iightcurves, caused by the different spacecraft 
orbits; we estimate that this is a ~ 10% effect. The remain- 
der is probably the result of calibration uncertainties, note 
that a 20% discrepancy between RXTE and ASCA has 
been reported by Yaqoob (2000). To overcome the prob- 
lem we fitted a simple power-law model to both ASCA 
and RXTE spectra, using only the 2-6 keV region of over- 
lap; we then adjusted the exposure times of the spectra 
until the fitted normalisations matched. In the following 
modelling we used the data ranges 0.8-9 keV (SIS), 0.8-11 
keV (GIS) and 2-25 keV (PCA). 

3.1. The basic model 

The accretion column of an intermediate polar is expected 
to consist of ~ 30 keV plasma at the accretion shock, cool- 
ing to < 1 keV nearer the white dwarf, and becoming opti- 
cally thick at some point. The temperature and density as 
a function of height are approximated by the Aizu model 
(Aizu 1973). Ramsay (2000) has fitted the RXTE spec- 
trum of V1062 Tau using an Aizu distribution of tem- 
peratures in a paper aimed at deducing the white-dwarf 
masses in the class of intermediate polars. Here, we are 
investigating the spectral changes over spin phase, and so 
use a simpler, but adequate, parametrisation involving a 
two-component plasma with temperatures of typically 1 
and 30 keV. 

The multi-temperature spectrum is modified by ab- 
sorption and scattering by the infalling accretion flow. We 
found that we needed both simple absorption and partial- 
covering absorption to reproduce the spectrum. We also 
obtained a better fit by introducing an iron edge fixed at 
7 keV, as would arise from strong ionized absorption (fits 
without this edge typically had \l poorer by 0.1). Lastly, 
we added a narrow line fixed at 6.4 keV to model the iron 
fluorescence arising from reflection, though for simplicity 
we did not include a reflection continuum. 

When fitting GIS, SIS and PCA spectra simultane- 
ously in 4 phase bins, the above model gave fits with \l 
typically 1. 1-1.2. If the model was optimised for each indi- 
vidual spectrum the resulting \l dropped to a 1, suggest- 
ing that calibration uncertainties are largely responsible 
for the higher m fh e simultaneous fits. 

3.2. A variable-absorption model 

We present in Table 1 a model which varies the absorption 
to reproduce the spin pulse, fixing the temperatures and 
normalisations over spin phase. The fitted temperatures 
(0.16 and 35 keV) are ill-constrained and, since we are 
not fitting with a full column model, we treat them as 
indicative only. Also fixed over spin phase are a simple 


Table 1. The variable-absorption model. I he basic model is 
a 0.16-keV MEKAL (normalisation 0.74) plus a 35-keV MEKAL 
(normalisation 0.023), a simple absorbing column of 1.7 x 10 22 
cm'" 2 , a 7.0-keV edge of t = 0.12, and a 6.4-keV line. The 
columns and covering fractions (/) of the partial absorption are 
tabulated. The fit to the S1S+GIS+PCA data simultaneously 
in 4 phase bins gave \ 2 V = 1.16 (^ = 1548). 

Ma x Fall Min Rise 

nH (xlO 22 cm" 2 ) - 61 8.6 19 

f - 0.40 0.69 0.33 


nH (xlO 22 cm 2 ) - - >1° 3 



absorbing column of 1.7 x 10 22 cm** 2 , a 7.0-keV edge of 
optical depth 0.12, and a 6.4-keV line. 

From here we reproduce the spin pulse by changing the 
column and covering fraction of one partial-covering ab- 
sorber. Further, to model the spin-minimum bin we need 
an additional dense absorber with column >10 25 cm’ 2 
covering 0.13 of the source. Of course, given the high col- 
umn, this is equivalent to a 13% flux reduction. 

Thus we can model the phase-resolved ASCA A RXTE 
spectra using 7 phase-independent parameters, a further 2 
absorption parameters per phase bin, and the normalisa- 
tion reduction at spin minimum. The overall was 1-1®* 
Adding further parameters, such as allowing the simple 
column and the optical depth of the / .0-keV edge to vary 
over spin phase, reduces the Xu 1-09, but makes little 
difference to the model as a whole. 

Modelling the spin pulse with phase-variable absorp- 
tion can be interpreted in the standard accretion-curtain 
model for IPs (Rosen etal. 1998; Hellier etal. 1991; Kim 
k Beuermann etal. 1995; Hellier etal. 1996) in which the 
upper pole points towards the observer at pulse minimum. 
The accretion flow then obscures the line of sight to the 
accretion regions, causing the absorption. In this model a 
partial absorber can be interpreted as either partial cov^ 
ering of a spatially extended emitting region, or the time- 
average of varying absorption. The 13% flux reduction at 
spin minimum then implies that part of the curtain is 
dense enough to be opaque, and that this covers the emit- 
ting region at spin minimum, either partially or for a frac- 
tion of the time. The column of >10 25 cm’ 2 is consistent 
with looking down an accretion column end-on (see, e.g., 
the column estimates in Hellier etal. 1996). 

Alternatively, it is possible that the disappearance of 
the lower pole over the white-dwarf limb at spin minimum 
might contribute to the flux reduction. In the simplest 
model of a symmetric dipole and zero shock height, the 
appearance of one pole compensates for the disappearance 
of the other. However, an asymmetry or a non-zero shock 
height allows both poles to be visible when the upper pole 
points away from us (spin maximum), with only one pole 
visible when the upper pole points towards us (e.g., Mukai 
1999). 


Table 2. The variable-normalisation model. It consists of a 1.5- 
keV mekal plus a 34-keV MEKAL, a simple absorbing column 
of 1.5 x 10 cm , another 10.1 x 10~ 2 cm" 2 column covering 
d% of the source, a 7.0-keV edge of r = 0.18, and a 6.4-keV 
line. The normalisations of the two plasma components are 
tabulated. The fit gave \l = 1.24 (v = 1550). 

Max Fall Min Rise 

3'1-keV norm (xlO -2 ) 2.56 1.95 2.15 2.36 

1.5-keV norm (x!0~ 2 ) 3.47 1.39 0.0 1.86 


3.3. A variable-normalisation model 

As an alternative to varying the absorption, we have ex- 
plored whether the spin pulse can be modelled by chang- 
ing only the normalisation of the spectral components, as 
would occur if parts of the accretion region were passing 
over the white-dwarf limb. Even if the absorption is not 
phase- varying, a partial-covering absorber of ~ 10 23 cm -2 
is still required to model the spectrum: multi-temperature 
models with no absorption, or with much lower absorp- 
tion, give xt > 2.4. We thus use the same basic model 
as above, a two-temperature plasma with partial-covering 
absorbers, an iron edge, and an iron fluorescence line. 

From here, the spin pulse is reproduced by modulat- 
ing the low-temperature plasma by 100% while the high- 
temperature plasma is modulated by 30% (see Table 2). 
This model gives xt = 124, somewhat worse than the 1.16 
obtained with the variable-absorption model. 

We can ask, further, whether this model makes physi- 
cal sense. For instance, it is unlikely that a 10 23 cm -2 ab- 
sorber covering 71% of the source would not change over 
phase; such absorption can only come from the magneti- 
cally trapped accretion flow, and this would wheel round 
on the spin cycle, coming into and out of the line of sight. 
Second, producing a spin-pulse by changing the normalisa- 
tions requires that parts of the accretion regions be passing 
over the white-dwarf limb at all spin phases (see Table 2), 
and this requires large, extended accretion regions, com- 
parable in size to the white-dwarf surface. However, this 
conflicts with the fact that the low-temperature compo- 
nent is 100% modulated, since, even if there is some asym- 
metry between the upper and lower poles, it is inconsistent 
that the accretion regions are (1) comparable to the white 
dwarf in size, (2) disappear completely over the white- 
dwarf limb, and (3) are not replaced by the appearance of 
the opposite accretion region. Further, the requirement for 
large accretion regions conflicts with the finding that they 
cover only <0.002 of the white-dwarf surface in XY Ari 
(Hellier 1997). For all of the above reasons, we consider 
that reproducing the spin pulse by varying normalisations 
alone is not viable. 

Of course, we cannot rule out the possibility that oc- 
cultation plays a secondary role, since allowing both ab- 
sorption and occultation to vary over spin phase produces 
a wide range of acceptable but non-unique fits. 


4. Discussion 

To model the X-ray spin pulse of V1062 Tau we require 
dense partial absorbers varying over spin phase. At sp,n 
minimum the absorber is sufficiently dense as to be com- 
pletely opaque. These results are consistent with absorp- 
tion in the accretion curtains of material falling onto the 
magnetic poles. However, owing to ambiguity in the model 
fits, we cannot determine whether, in addition, parts of the 
accretion regions are passing over the white-dwarf limb. 
The need to average over substantial ranges of spin-phase, 
and the complexity of the models when both absorption 
and the normalisations of multi-temperature plasmas are 
allowed to change, means that unique solutions are not 
obtained. In reality, we should also have included a reflec- 
tion continuum, adding further phase- varying parameters 
ancl making the ambiguity worse. 

Thus, even with combined ASCA and RXTE data, 
the complexities of apparently simple sinusoidal pulsations 
are such that spectral analysis is a poor diagnostic of the 
accretion flow, beyond the statement that dense (>10 23 
cm ), phase- variable absorption is a primary factor. 
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